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Abstract 

Using extensive computer simulations, the behavior of the structural modes — more precisely, the 
eigenmodes of a phantom Rouse polymer — are characterized for a polymer in the three-dimensional 
repton model, and are used to study the polymer's dynamics at time scales well before the tube 
renewal. Although these modes are not the eigenmodes for a polymer in the repton model, we show 
that numerically the modes maintain a high degree of statistical independence. The correlations in 
the mode amplitudes decay exponentially with (p/N) 2 A(t), in which p is the mode number, N is 
the polymer length and A(t) is a single function shared by all modes. In time, the quantity A(t) 
causes an exponential decay for the mode amplitude correlation functions for times < 1; a stretched 
exponential with an exponent 1/2 between times 1 and tr ~ A^ 2 , the time-scale for diffusion of 
tagged reptons along the contour of the polymer; and again an exponential decay for times t > tr. 
Having assumed statistical independence and the validity of a single function A{t) for all modes, 
we compute the temporal behavior of three structural quantities: the vectorial distance between 
the positions of the middle monomer and the center-of-mass, the end-to-end vector, and the vector 
connecting two nearby reptons around the middle of the polymer. Furthermore, we study the 
mean-squared displacement of the center-of-mass and the middle repton, and their relation with 
the temporal behavior of the modes. 

PACS numbers: 36.20.-r,64.70.km,82.35.Lr 



2 



I. INTRODUCTION: POLYMER REPTATION AND THE REPTON MODEL 



The dynamics of a coarse-grained polymer has two extremes in polymer physics. One 
extreme features the phantom Rouse model jl, 2|, wherein the polymer is described as a chain 
of beads serially connected by harmonic springs. A polymer in the phantom Rouse model 
does not feel its surroundings, and not even itself. Comprised of iV + 1 beads with spatial 
locations r*j (i = 0, 1, ... , N), a polymer of length N in the Rouse model lends itself to a 
complete analysis of its dynamics in terms of the eigenmodes indexed by p (p = 0, 1, . . . , N), 
correspondingly called the Rouse modes, whose amplitudes are given by 

1 - 

x pW = WTi ^ cos [ap{l + 1/2)] m ' (1) 

8=0 

in which a p = pir/(N + 1). The other extreme features Rubinstein's repton model 
which describes the dynamics of a highly restricted polymer, in which the only dynamical 
freedom is a mechanism known as reptation , a snake-like motion in which stored 

length travels back and forth along the polymer, and can be created and annealed only at 
the polymer's ends. In the repton model, the polymer is composed of mobile points, called 
reptons, residing in the cells of a lattice, connected by bonds. A configuration of the polymer 
in the repton model in two dimensions is shown in Fig. [TJ The neighboring reptons can 
reside either in the same cell or in the adjacent cells. The configuration of the polymer's 
contour in space, defined by the serially connected cells occupied by the reptons, is that 
of a random walk, and is not allowed to undergo any change, other than by the motion of 
one of its ends. Each repton, except the end ones, is therefore constrained to move along 
the contour of the polymer. The repton model is extensively discussed and illustrated in 
the literature Q B. Duke has extended the repton model to study gel electrophoresis Q, 
by adding an electric field which acts on the charged reptons; also this model has been 



studied extensively |6|, [8, 9]. A highly efficient simulation approach of the repton model was 
developed based on bit-operations [8]; the present work uses this approach as well. 

In this paper we characterize the analogs of the Rouse modes for a polymer consisted of 
N + 1 reptons in the repton model in three spatial dimensions, having denoted the spatial 
locations of the reptons by r*j (i = 0, 1, . . . N), as in Eq. ([IT) . Henceforth, these are simply 
referred to as the modes. They are not the eigenmodes of the polymer in the repton model. 
We express several dynamical variables for the polymer in terms of the modes. We verify 



the time-correlations of the modes and these dynamical variables by high-precision computer 
simulations. 

In the simulations the repton positions are the natural variables. For theoretical analysis 
however, the distance between successive reptons, or the "link variables", defined as 



yi = n- n-i 



for % — 0, 



N. 



(2) 



are a more convenient choice. By construction, the link variables are either a unit vector, 
or zero. In the former case, the link is taut, while in the latter case the link is slack and 
it represents an element of "stored length". The motion of the polymer can be viewed 
as diffusion of stored length along the chain (instead of the movements of the individual 
reptons). 




FIG. 1: A two-dimensional lattice representation of the repton model. The polymer is represented 
by points (reptons) that reside in the cells of the lattice, and are connected by bonds. The rules 
governing this model are described in the text. 

Few exact results are known for the repton model, all concerning the asymptotic time 
regime, governed by the overall diffusion constant D of the polymer. It was already proposed 
yy Rubinstein j^], that D scales as iV -2 . Later it was found by van Leeuwen and Kooiman 



lOfl that in the repton model on a rf-dimensional hypercubic lattice, D = l/[(2d + 1)N 2 ]. 
This result was shown to be exact by Prahofer and Spohn {ll| and by Widom and Al- 
Lehyani [12]. In this paper, we study the repton model on a cubic lattice. For reasons of 



computational efficiency, the rates for extraction and retraction of the ends are modified 
to obtain a dimensionality-independent density of stored length, tuned to that of the one- 
dimensional model, so that D = 1/(3N 2 ). We are specifically interested in the polymer's 
dynamics over intermediate time regimes, for which, so far, only simulations can reveal the 
behavior. 

Before the diffusive regime is reached, the repton model has various distinct time regimes. 
First, at very short times t < 1, the dynamics is governed by individual hops of reptons, 
which are spatially separated. Next, there is one intermediate regime 1 < t < tr ~ N 2 , 
where tr is the time at which the fluctuations in the density of stored lengths decay. This 
is followed by a second intermediate regime tr < t < ~ N 3 , where is the well-known 
tube-renewal time after which the initial information on the chain has been forgotten. 

Our basic quantity of interest is the equilibrium correlation function 



in the above time regimes, where the angular brackets denote an average over the equilibrium 
ensemble of polymer conformations. More specifically, we put the middle monomer at t = 
at the origin. Then a polymer configuration can be build up from independent link variables 
i/i. The link may be taut or slack, with probabilities 2/3 and 1/3. A taut link variable may 
assume one of the 6 possibilities corresponding to the six directions on the cubic lattice, 
each with equal weight. 

Since the modes are not eigenmodes of the polymer in the repton model, they do not 
decay exponentially in time. In Sec. [Til we show that the cross-correlation functions for p ^ q 
vanish for t — 0, and remain small at later times. In Sec. II III we study the autocorrelation 
functions C pp (t). These autocorrelations are found to decay exponentially at short times 
t < 1 and at long times t > tr, but in the intermediate regime 1 < t < tr they decay 
as stretched exponentials in time, with exponent 1/2. In all these regimes, the exponent 
has a power-law dependence on mode number p and polymer length N. In Sec. IIV| using 
the mode amplitude autocorrelation functions C pp (t) and neglecting cross-correlations, we 
reconstruct other time-dependent correlation functions, in particular those of the vector 
between the positions of the middle repton and the center-of-mass, the end-to-end vector, 
and a spatial vector connecting two nearby reptons around the middle repton. We compare 
these analytic results with direct computer simulations. In Sec. [V] we evaluate the mean- 




(3) 
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squared displacement of the center-of-mass and of the middle repton in the above time 
regimes, and compare them with direct computer simulations. We end the paper with a 
discussion in Sec. [VD 

II. BEHAVIOR OF CORRELATIONS BETWEEN DIFFERENT MODES 

For phantom chains in the Rouse model, the modes as defined in Eq. (CQ) are strictly 
independent and orthogonal, i.e. the correlations C pq (t) as defined in Eq. are zero for all 
times, if p ^ q. This property is part of the motivation why in the study of the dynamics 
of a polymer, it is very convenient to express the quantities of interest in amplitudes of 
the structural modes. In this section we will show analytically that at t = 0, also in the 
repton model the same structural modes are orthogonal. We will also show that correlations 
between even modes p and odd modes q are strictly zero at all times. Next, with computer 
simulations, we will show that for time t > cross-correlations between the modes p ^ q, in 
which p and q are both even or both odd, start to develop. 

A. Behavior of C pq (0) 

First, we turn to the correlation functions C pq (t) for time t = 0, which only require the 
evaluation of equilibrium averages. For this purpose it is convenient to provide the relation 
between the mode amplitudes and the link variables which have the simple equilibrium 

averages 

(ViVk) =I(p/d)5 jk , (4) 

where I is the unity tensor, and p is the equilibrium density of the taut links, which equals 
p = 2/3 in our simulations. 

We express the position of the i-th repton w.r.t. that of the middle repton [for 
simplicity we take N even all throughout this paper; if one is interested in odd N polymers, 
one can take as a reference point r m = | (r(N-i)/2 + r(jv+i)/2j , which has the same symmetry] 
and the link variables as 

3 N/2 

rj = f N/2 + Yl for > N / 2 )i ?i = ?n/2 ~ Vi for (j < N/2). (5) 

i=N/2+l i=j+l 
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In the summations to come we will use the following relation for p^O 

N 



^cos[a p (j + l/2)] 

j=k 



sm(a p k) 

2.Sr) 



(6) 



with s p = sin(a p /2). Note that Eq. ([6]) yields for k — 0. This implies that the middle 
repton location r/v/2 in Eq. (j3J) does not contribute in the mode amplitude expression for 
p 7^ 0. Thus, using Eq. d^J), and further, upon interchanging the summations over j and k, 
for p 0, we have 



1 - 

+ I s, 



sm(a p k) y k . 



(7) 



2(iV , ^ p ^ 

Equation ([7j) cannot be applied to the center-of-mass mode p = as sq = 0. For p = 0, with 
the help of Eq. (jSJ), the corresponding formula becomes 



X = fzv/2 + 



iV + 1 



AT/2 

XI i (yx+i-j -Vj) 



(8) 



While the mode amplitudes for p 7^ are expressed solely in terms of the link (or structural) 
variables, the center-of-mass mode needs an additional (translational) coordinate rjv/ 2 - 
Equation (J7]) can be inverted using the relation 



N 



^ sin(a p £;) sin(a q k) = — y 1 5 OT 
fc=i 



and its inverse 



TV 



^ sin(a p A;) sin(a p /) = 1 S M . 



(9) 



(10) 



P =i 



These results lead us to the inverse relation for the link variables in terms of the mode 
amplitudes 

TV 

" ~ (11) 



Vk 



-4 s p sm(a p k) X p . 



Equation flllj) reveals that the mode amplitudes for p 7^ and the link variables are fully 
equivalent since they can be transformed into each other. 

From the expression (JHJ) and the average (j3J) it follows that the mode amplitudes for p 7^ 
and q 7^ are orthogonal for t = with the inner product 



C Pg (0) 



P 



l(N + l)sl 5pg - 



(12) 
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B. Behavior of C pq (t) for p, q > and t > 

In Fig. [2] we plot the normalized cross-correlation 

n Pq {t) = c Pq {t)/^c pp {t)c qq {t) (is) 

for a number of p and q and several times. In the upper left panel one observes the absence 
of cross-correlations for t = 0, as derived in Eq. (1121) . The next clear observation in Fig. [2] 
is the absence of cross-correlations between the even and odd modes. This is a consequence 
of the inversion symmetry of the polymer chain: inversion is a renumbering of the reptons 
from N to and the inversion of the link variables, or notationally, 

Vi <-» -m-i+i- ( 14 ) 

The amplitudes of the even modes are invariant under this operation, while the odd mode 
amplitudes change sign. Since the dynamical evolution commutes with the inversion op- 
eration, the symmetry of the even mode amplitudes and the change of sign for the odd 
mode amplitudes are preserved under time evolution. In the lower-right panel the inversion 
symmetry appears to be violated for larger p and q, but this is actually a consequence of 
the numerical inaccuracy of the simulation data, which is magnified once the normalizations 
involved in rj pq become very small. 

Figure [2] also reveals that for t > 0, C pq (t) is no longer strictly proportional to 5 pq . 
Nevertheless, the off-diagonal correlation functions are small in comparison to the diagonal 
ones. 



III. BEHAVIOR OF C pp (t) 

For the properties for C pp (t) we can anticipate the following behavior. 

(i) At short times t < 1, reptons typically make zero or one move, and the movements 
of individual reptons are uncorrelated. In an earlier work on lattice polymers with 
exactly the same equilibrium properties as in the repton model, but less restricted 
dynamics [l^|, we found that the decay of the structural modes in this short-time 
regime was characterized by C pp (t) = C pp (0) exp [—AipH/N 2 ], for some constant A\. 
Since the moves allowed in the repton model are a finite fraction of the moves in this 
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FIG. 2: Checkerboard plot for rf pq (t) for N = 400 at four different times; upper left: t = 0, 
upper rig lit: t = 8,000, lower left: t = 64,000, and lower right t = 10 6 . The data are manually 
symmetrized around the diagonal. The grayscale used is shown in the scale bar. 

earlier work, we expect the same behavior to hold, be it with a modified value of the 
parameter A 1 . 

(ii) In the regime t > 1 reptons have the time to attempt multiple moves. Right after a 
repton has moved in the "upstream" direction, the reverse "downstream" move is guar- 
anteed to be possible, while another "upstream" move requires the presence of another 
slack link upstream, and therefore has a significant probability to be ruled out. Hence, 
the nature of the dynamics changes at t ~ 1. From a curvilinear perspective, the 
motion changes from ordinary diffusion to "single-file diffusion", in which the squared 
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FIG. 3: Behavior of C pp (t) for t > 0. (a) for t < 1, C pp (t) = C pp (0)exp[-A lP 2 t/N 2 } for some 
constant A 1 of magnitude O(l). (b) for 1 < i < tr ~ iV 2 , C pp (i) = C pp (0) exp[-A 2 pV /2 /iV 2 ], 
while for i > tr, C pp (t) = C pp (0) exp[— A 3 p 2 t/N 3 ], with and ^3 also of magnitude 0(1). 

curvilinear displacement increases ~ \rt. Continuity of the curves near t ~ 1 then 
leads to an expected behavior C pp (t) = C pp (0) exp [— A 2 p 2 \/i/N 2 ~\ . 

(iii) Slack links can move freely up and down the chain, in unit steps, without con- 
straints. This diffusive behavior along the chain with a curvilinear diffusion constant 
of ~ 1, combined with the polymer length ~ N, yields a characteristic time scale 
tr ~ iV 2 beyond which the distribution of slack links becomes statistically indepen- 
dent. Without such correlations, the curvilinear dynamics ceases to be anomalous and 
becomes ordinary diffusion once more. Continuity of all curves around tr then yields 
C pp (t) = C pp (0) exp [-A 3 p 2 t/N 3 ]. The data reveal that A u A 2 and A 3 are all of 0(1). 

Having combined (i— iii) , we present the simulation data for several p and N values in Fig. 
[3j The data collapse confirms that for t < 1 

C pp {t) = C pp (0) exp[-A lP 2 t/N 2 } (15) 

for some constant A\ of magnitude 0(1). Further, for 1 < t < tr ~ iV 2 

C pp (t) = O pp (0) exp[-A 2 p 2 t^ 2 /N% (16) 

while for t > tr 

C„(t) = O pp (0) exp[-A 3 p 2 t/N 3 }. (17) 
10 



Having combined Eqs. (I15til7p . we see that C pp (t) is expressed in general terms as 
C pp (0) exp[—(p/N) 2 A(t)} with a time-dependent quantity A(t); which we will encounter in 
the coming sections. At short and long times, this function A(t) is linear in time, with at 
short times a length- independent prefactor A\ and at long times a prefactor A 3 /N, inversely 
proportional to polymer length; in the intermediate regime 1 < t < Tr ~ iV 2 , the function 
A(t) increases proportional to with a length- independent prefactor A^. Note that at all 
times A(t) is independent of the mode number p. 



IV. THREE STRUCTURAL QUANTITIES IN TERMS OF MODE AMPLITUDES 

All structural quantities which can be expressed in terms of the link variables, are also 
expressible in terms of the mode amplitudes. In the previous section, we obtained an ap- 
proximate description of the behavior in time of the mode amplitudes: we neglect cross- 
correlations at all times, and assume that a single function A(t) determines the autocorre- 
lation of all mode amplitudes. In this section, we will confront this approximate description 
with direct simulations. For this purpose, we consider three such structural variables: the 
vector So between the center-of-mass and the middle repton, the polymer's end-to-end vector 
Si, and the vector S2 connecting the two reptons next to the middle repton. 



A. Structural variables So, Si and S2 expressed in terms of the modes amplitudes 



The distance So between the center-of-mass and the middle repton, defined in Eq. (jH 
is given by 

1 



So — Ro — r N/2 



N + l 



N/2 

im+i-j-yj) 
3=1 



(18) 



in terms of the link variables. Its mean-squared equilibrium-average follows from Eq. (j^J) as 



(So 2 ) 



2X P N(N + 2) 



(19) 



12(iV+ 1) ' 

Having expressed the link variables in terms of the mode amplitudes further yields the 
expression 

N 

(20) 



So = 2 J2 (-!) P/ % 

p=2,4,..- 
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Similarly, the end-to-end vector Si of the polymer is expressed as 

TV 

Si = r N - f = ^ sT fc . (21) 
fc=i 

in terms of the link variables. Its mean-squared equilibrium-average is given by 

(S?) = P N, (22) 
and the expression of S± in terms of the mode amplitudes takes the form 

JV 

5i = -4 ^ c p X p . (23) 

p=l,3,- 

with c p = cos(a p /2). 

Note that the end-to-end vector is a special case of the vector between two arbitrary 
reptons m and n 

m 

£ ^ (24) 



fc=n+l 



For a characteristic case, we consider the vector S2 between the two reptons next to the 
middle repton as our third structural variable, i.e., 

S 2 = r N /2 + i - rjv/2-1, (25) 

with the mean-squared equilibrium-average 

(S 2 2 > = 2p; (26) 

5*2 is related to the mode amplitudes as 

N 

f N/2+1 ~ r N /2-i = -8 £ sin(p7r/2) c p s p X p . (27) 

p=i 

We summarize the expressions of So, ^1 an d S2 [Eqs. (120|) . (123]) and (127|) ] in the form 

v 

Sj = S p ^X p . (28) 
P =i 

The equilibrium averages can also be evaluated via the averages of the mode amplitudes, as 
provided in Eq. ffT2l) . This leads one to the evaluation of sums of the type 

N 

^ 2 ) = WtI)E(^ p ) 2 , (29) 
' p=i 
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which can be found in Ref. 13] . 



Although the mode amplitudes for p ^ are equivalent to the link variables, one must 
realize that the domains of possible values are quite different. The domain of possible jjj 
consists of 2d + 1 = 7 values for each j in three dimensions (d — 3), whereas the domain of 
possible values of the X p is complicated. Each X p assumes (2d + 1)^ values (multiplicities 
included) and the possibilities are inter-related. 

B. Dynamical properties of the structural variables S 

Whereas the equilibrium averages of structural properties are easily evaluated using the 
link variables, their time- dependent correlators have to be evaluated through the use of 
temporal behavior of the correlations in the mode amplitudes. 

Using the general form Eq. (I28j) . we find, for the correlators, that 

N 

(Si(t) ■ Si(0)) = S P ,iS q ,iC M (t). (30) 

p,q=l 

For the time dependent C pq (t) we use the approximation discussed in the previous section, 
leaving out the cross-correlations and employing for the diagonal terms the generic form 

C pp (t) = C pp (0)exp[-(p/N) 2 A(t)] = 8{N P +1)s2 exp[-(p/N) 2 A(t)]. (31) 

In other words, 

(Si(t) .£(0)> = j^-rJ2( S P^ s p) 2e M-(p/N) 2 A(t)}. (32) 
+ p=i 

We note that the t = values are given as sums over the ratios (S p ^/s p ) 2 , which are most 
easily evaluated through their relation with the link variable representation. 

For the finite chains that we consider here, numerical evaluation of the sums presents no 
problems. The behavior can be generally characterized by three regimes. 

1. The initial regime, in which A(t) remains of order unity and all terms in the sum over 
p participate. In this regime the deviation from the initial value is proportional to 
A(t). 
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2. The intermediate regime from 1 < t < Tr where A(t) grows as \/i with time. Then 
the small p modes contribute to the sum, but the exponential effectively cuts the sum 
off at a value of p ~ N/ a/ ' A{t). 

3. The asymptotic regime, starting from the Rouse time tr. Then A(t) is of order 
1/N and the sum is effectively restricted to its first term with exponentially small 
corrections for the other p. 

In order to check the approximation fl32l) we performed computer simulations for chains 
with length iV = 400, and stored at regular times the instantaneous values for the struc- 
tural variables So, Si and S2. These observables show a dynamical behavior with variations 
spanning over many decades. For practical purposes, we therefore performed independent 
simulations with n^t = 1, 10 2 , 10 4 , 10 6 or 10 7 elementary moves between successive measure- 
ments, corresponding to time steps of At = n& t / (iV+1). For each time step, we generated 64 
sets of 10 5 measurements. The results of these independent simulations were then combined 
into a single curve in Fig. [5] (see also the next paragraph), which occasionally shows little 
"hiccups" at the times where the results of independent simulations are joined together. 




FIG. 4: The correlation functions Si (t) = 0i(t) ■ Si(0))/(Sf) for i = 0, 1 and 2. 

In Fig. H] we plot the measured value of the correlators for iV = 400 and the calculated 
values from Eq. (132|) . The agreement is very good for i = and i = 2; for i = 1 the calculated 
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FIG. 5: Behavior of the mean-squared displacement of the center-of-mass, (Ar 2 m (t)). (a) for t < 1, 
<Ar 2 n (i)> ~ i/JV. (b) for 1< t < ~ iV 2 , <Ar c 2 m (t)> - Vt/iV, while for t > t r , <Ar c 2 m (t)> ~ t/N 2 . 

values are about a factor of two smaller than the measured values in the intermediate range. 
Given the very wide range in time and values of these correlators the correspondence is still 
impressive. Since the main approximation is the neglection of the cross-correlations in the 
structural modes, this is indirect proof that although the number of cross-correlations is 
very large, the main contribution comes from the autocorrelations. 



V. MEAN-SQUARED DISPLACEMENT OF CENTER-OF-MASS AND MIDDLE 
REPTON 

Besides the dynamic behavior of structural properties of the polymer, also the dynamical 
behavior of its position is of interest. Two key observables that characterize this, are the 
mean-squared displacement of the center-of-mass, and of the middle repton. 

We performed computer simulations of polymers with length N = 100, 200, 400 and 800, 
and stored at regular times the positions of the middle repton and the center-of-mass. For 
each length, we performed independent simulations in which the time between updates 
of the positions was n At = 1, 10 2 , 10 4 , 10 6 or 10 7 elementary moves, corresponding to time 
differences of At = n At /N. For each length and time difference, we performed 64 simulations 
of 10 5 time steps each. 

We first show the results for the mean-squared displacement of the center-of-mass and 
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FIG. 6: Behavior of the mean-squared displacement of the middle repton, (Ar^y 2 (t)). (a) for 
t < 1, (Ar 2 N/2 (t)) ~ t. (b) for 1 < t < t r ~ iV 2 , (Ar^ /a (t)) ~ t 1 / 4 . (c) for r fi < t < r Cj 



(Ar^ /2 (i)) ~ ^/pf, and for t > r c , (Ar^t)) ~ t/iV 2 . 

the middle repton of the polymer, and then we discuss how the two are related. 

The computer simulations show that the mean-squared displacement of the center-of-mass 
of the polymer behaves as follows: 



t/N 



t < 1 



(IMt) - Mo)] 2 ) = (ALW) ~ 



t 1/2 /N l<t<r R ~N 2 



(33) 



t/N 2 



t>TR 



The data corresponding to Eq. ( l33l) are shown in Fig. |5j 

For the mean-squared displacement of the middle repton we find from the simulations 



([r N/2 (t) -r N/2 (0)} 2 } = (A 2 N/2 (t)} 



l<t<r R ~N 2 

TR<l<T d ~ N 3 
t > T d 



(34) 



t 

t l/4 
y/t/N 

[t/N 2 

The data corresponding to Eq. (134p are shown in Fig. |6j Note that the microscopic formu- 
lation for the mean squared displacement of the middle monomer is that of a Generalized 
Langevin Equation |l5| . 

The two mean-squared displacements are related in the following manner 



(AL(*)> - = 2 <^W ■ - 3>(0)]> + - 5 (0)] 2 ), 



(35) 
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since So is the vectorial distance between the center-of-mass and the middle repton. The 
last term in Eq. (135|) can be related to a structural correlation function using the relation 

(S 2 (t)) = (S 2 (0)}. (36) 

This equality holds because the structural distribution is not affected by the position of 
the middle repton at t — 0. So we may evaluate this terms as 

([S (t) - S (0)f) = 2[<So 2 > - (S (t) ■ S (0))]. (37) 

The correlator of So has been worked out in Sec. IIVBI and the equilibrium value in Eq. 

din]). 

The first term (1351) . being a cross-correlation between a translational and a structural 
variable, cannot be computed using the structural modes only. It is interesting to remark 
that 

(r N /2(t) ■ [S (t) + So(0)]) = 0, (38) 

since the sum So(t) + £o(0) is time-reversal invariant, while f^^it) changes sign under 
time reversal. We have verified that the simulations obey this consequence of time-reversal 
invariance. We therefore have 

(r N /2(t) ■ [S (t) - S Q (0)]) = -2(f N/2 (t) • 5 (0)). (39) 

It is clear that the right hand side of Eq. fl39l) is negative because the middle repton has 
a tendency to move in the direction of the instantaneous position of the center-of-mass. It 
is worth noting that this tendency has a persistent effect on the value of the correlation 
function. Similarly, it can be shown that the center-of-mass also has a tendency to move 
towards the middle repton! This means that a reptating polymer experiences more than 
only internal forces — indeed, part of the forces a reptating polymer experiences originate 
from the dynamical constraints that prevent the polymer from moving sideways; giving rise 
to the phenomenon that mode p = is not independent from the other non-zero modes. In 
Fig. [7] we show simulation results for the first term in the left-hand side of Eq. ( I35jl . Note 
that the correlation function saturates asymptotically to a value close to N/9. 

Let us now remark here that the t 1//4 behavior of the middle repton is a consequence of 
two separate aspects of the middle repton's motion: (i) that the middle repton simply moves 
along the backbone of the polymer (curvilinear motion). For t < 1 the middle repton moves 
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FIG. 7: The cross-correlation function q(t) = — (fjv/2(*) " [<Sb(£) — •S'o(O)])- 

diffusively along the backbone of the polymer. For 1 < t < Tr ~ A^ 2 its motion corresponds 
to that of single-file diffusion, which yields the curvilinear mean-squared displacement of 
the middle repton (Au 2 N ^ 2 (t)) c , where the subscript 'c' denotes an average for a given initial 
polymer configuration, increasing in time as t 1 / 2 . For t r < t < r d ~ iV 3 the motion of the 
middle repton becomes diffusive along the backbone, and its curvilinear motion becomes 
diffusive in time, i.e., increases ~ t/N, and beyond t r the backbone of the polymer is 
completely renewed in physical space, (ii) The backbone configuration of the polymer, 
from one end to the other, is that of a random walk in physical space. Hence, the real- 
space mean-squared displacement of the middle repton (Ar 2 V//2 (t)) averaged over all polymer 
configurations equals ' (Au 2 N ^ 2 (t)) c ; i.e., (Ar 2 v ^ 2 (t)) ~ t 1 / 4 for 1 < t < tr, and then ~ \/t/N 
for r R <t< T d . 

Given the above discussion, one can relate the mean-squared displacements of the middle 
repton and the center-of-mass in the following manner. First, we ignore the end-effects, i.e., 
assume that the motion of all the reptons progress along the polymer's backbone, which is 
fixed in space. Secondly, for a given backbone of the polymer fixed in space we imagine 
an interior repton i of the polymer at t — 0. After time t, the repton would have moved 
to a different location on the backbone, which was occupied by another repton j at t = 0: 
clearly, j is a function of i and t, i.e., j = ji(t). In fact, the real-space displacement Afj(i) 



18 



for repton % is then written as, using Eq. (j2J), 

kit) 

Af l (t) = J2m(o)- (40) 

k=i 

Thereafter, in terms of the link variables {yk} s , the real-space displacement of the center- 
of-mass for this polymer in time t is given by 

1 N+l , N+lji(t) 

Af -W = (jvTT) £ AfiW = (jvTT) £ £ «<»■ («) 

Thirdly, as we observe that for 1 < t < tr, (Ar 2 (i)) c ~ i 1 / 2 , it becomes clear that 
— i| ~ t 1//4 , meaning that in Eq. (14 ip each distinct yk(0) appears ~ t 1 / 4 times. With 
(yfci(0) ■ yfc 2 (0)) = <5fcifc 2 5 this implies that (r 2 m (t)) ~ y/i/N. Note that this argument can 
be used again to show that (r 2 m (t)) ~ t/N for t < 1 and also that (r 2 m (t)) ~ t/A^ 2 for 

< t < r^. Indeed, as for the pure time behavior is concerned, this argument shows that 
for a general reptation motion of a polymer, if the real-space mean-squared displacement 
of the middle repton increases in time as t a for some a, then the real-space mean-squared 
displacement of the center-of-mass has to increase in time as t 2a /N. 

VI. DISCUSSION 

By extensive simulations we have established scaling properties of the structural modes 
for a polymer in the repton model. It turns out that their autocorrelation functions can be 
written as a function of the scaling variable p/N, where p is the mode number and TV the 
length of the chain, and a function A(t) which is shared by all modes. This latter function 
causes exponential decay at short times t < 1 and long times t > Tr ~ iV 2 , and causes an 
intermediate regime of stretched exponential decay. 

Although the notion of structural modes stems from a freely moving polymer, held to- 
gether by harmonic forces, it also turns out to be useful for the dynamics in the repton 
model, for which the chain is restricted to move along its contour. The reasons are, apart 
from the above mentioned scaling properties, the following. 

1. The structural modes form a complete basis for the structural quantities, i.e. any 
linear function of the link variables can be written as a superposition of structural 
modes. 
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2. The equal-time correlations between the structural modes are strictly orthogonal and 
they preserve this orthogonality to a high degree for correlations at different times. So 
for practical purposes cross-correlations between structural modes may be ignored. 

3. This orthogonality allows to write the correlations in time between structural quanti- 
ties as a sum over the autocorrelation functions of the structural modes. 

4. The motion of the center-of-mass is correlated with that of the middle monomer. This 
stems from the fact that a reptating polymer experiences more than only internal 
forces — part of the forces a reptating polymer experiences originate from the dynam- 
ical constraints that prevent the polymer from moving sideways. The effect of this 
correlation is captured in the phenomenon that mode p = is not independent from 
the other non-zero modes. 

The mean-squared displacement of the center-of-mass shows also changes in its behavior 
at the same crossover times t ~ 1 and t = tr ~ iV 2 : from an initially diffusive regime 
in which it grows as t/N similar to a free Rouse chain, via an intermediate regime with 
anomalous diffusion in which it grows as yJt/N, to an asymptotic regime with reentrant 
diffusive behavior where it grows as t/N 2 . 

The asymptotic regime of the mean-squared displacement of the middle repton equals 
that of the center-of-mass, as it should to preserve integrity of the chain, but this regime 
is only entered after the tube renewal time ~ N 3 . It has an initially diffusive regime 
in which it grows linearly with time, also similar to the free Rouse chain. In contrast to 
the center-of-mass, it shows however two intermediate regimes in which it grows as t 1//4 and 
\Jt/N, respectively. We have presented arguments that for general reptation motion of a 
polymer, a real-space mean-squared displacement of the middle repton scaling as t a for some 
a must be accompanied by a mean-squared displacement of the center of mass increasing in 
time as t 2a /N. 

Finally we note that the most interesting and revealing aspect of the mode analysis 
presented in this paper is the interplay between the translational degree of freedom and 
the the structural degrees of freedom, i.e., the cross-correlation between the position of the 
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middle repton and the vector connecting the center-of-mass and the middle repton. 
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